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Abstract 

We have examined theoretically the spontaneous thinning process of tip- 
suspended nanowires, and subsequently studied the structure and stability 
of the monatomic gold wires recently observed by Transmission Electron 
Microscopy (TEM). The methods used include thermodynamics, classical 
many-body force simulations, Local Density (LDA) and Generalized Gra- 
dient (GGA) electronic structure calculations as well as ab-initio simulations 
including the two tips. The wire thinning is well explained in terms of a ther- 
modynamic tip suction driving migration of surface atoms from the wire to 
the tips. For the same reason the monatomic wire becomes progressively 
stretched. Surprisingly, however, all calculations so far indicate that the 
stretched monatomic gold wire should be unstable against breaking, contrary 
to the apparent experimental stability. The possible reasons for the observed 
stability are discussed. 
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In the last few years the properties of atomic size metallic contacts have attracted great 
attention both for fundamental and for practical reasons. Experimental formation of these 
contacts has been achieved using different techniques and in diverse conditions ranging from 
room temperature under atmospheric pressure [l[], to low temperatures in ultrahigh vacuum 
0. The relationship between mechanical structure and electrical properties has been studied 
by simultaneously measuring load force and conductance ||. On the theoretical side, the 
mechanical structure and evolution of a tip-surface mechanical contact has been modelled 
using classical molecular dynamics (MD) simulations with many body potentials [§,0,0 , and 
also by ab-initio MD simulations. ||. The conductance of these necks has been calculated 
using free-electron , tight binding []8| and ab-intio schemes P,[T0[] . The relationship between 
mechanical load and conductance has been studied using different models like free electron 
nj, tight binding and LDA [|TJ . 



Very long tip-suspended nanowires of surprising regularity and stability were found in 
simulations close to the melting point and subsequently confirmed experimentally |T3 



Questions about the structural and thermal behavior of ideally infinite nanowires have re- 
cently stirred an upsurge of theoretical [IB,|IT| and experimental interest Jl9|. New observed 
structural features of long nanowires include a striking regularity ||19|| . Both crystalline and 
noncrystalline morphologies have been predicted [T7| and observed |19|| . Very recently the 
exceptional ductility and diffusivity of metallic gold has been put to use in creating the 
ultimately thin nanowire, consisting of just a single atomic strand [^,^,[1^]. Direct TEM 
imaging, |2(J in particular, has shown that a thin multi-shell wire suspended between two 



(110) tips will gradually and spontaneously thin down to smaller and smaller radii, until 
the two tips are eventually connected by just a 4-atom long monatomic wire. Two of these 
four Au atoms stretch out of a well-tapered tips, the remaining two are freely suspended 
in vacuum in an astonishingly stretched configuration. The stretching magnitude increases 
with time, the average atomic distances growing from about 3 to about 4 A, yet the wire 
is not seen to break for many seconds. In different measurements |2T],[B|, the wire was not 



directly imaged, but conductance plateaus close to and below G = 2e 2 /h indicated that 
the tip-suspended wire must have been monatomic and conducting. For that case, however, 
no structural information is available. These findings raise several questions, among them: 

i) why are the wires thinning down with time, and what is the origin of the stretching 
force which appears to act on the monatomic wire even in the absence of any mechanical 
tip withdrawal? 

ii) how can we understand the structural stability of such an exceptionally stretched 
monatomic wire, particularly in terms of its electronic structure? 

Hi) how can we explain the wire conductance? 

While we have not yet addressed Hi, the scope of this note is to point out that thermody- 
namic reasoning readily explains point i, whereas preliminary calculations and simulations 
performed with state-of-the-art methods totally fail to account for the wire stability issue 
of point ii. This is all the more surprising since these methods, tuned on the properties 
of atomic Au, are known to work quite well for bulk and surfaces of metallic Au ||22|| . In 
the following we will address first the thinning. Secondly we will describe our attempts 
at understanding cohesion and breaking of the infinite wire and of the tip-suspended short 
wire, and finally speculate on what might be the explanation for the unusual stability. 
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I. WHY DO SUSPENDED WIRES THIN DOWN, AND WHY IS THE 
MONATOMIC WIRE STRETCHED? 



All wires suspended between two tips appear to thin down steadily with time, on a scale 



which is visible within seconds. [|19| , p0|| . The rate of decrease moreover appears somewhat 
faster for thinner wires. Once the wire is monatomic, it becomes increasingly stretched with 
time, even though the tips are fixed relative to one another. 

We propose the following simple thermodynamical theory. For a sufficiently thick wire 
of radius R, the free energy per atom is: 

e(R) =e + 2 1 /Rp (1) 

where Eq is the bulk free energy (essentially identical to the cohesive energy), 7 is the surface 
free energy per unit area and p is the bulk density. Taking one atom from the wire to the 
tip, the free energy decreases, and in turn this represents an effective suction force exerted 
by the tips on the wire atoms. In presence of surface diffusion, the wire surface atoms 
will drift towards the tips, and the wire radius will thus decrease. For a surface diffusion 
coefficient D s and an Einstein-related surface mobility p s = D s /ksT, an atom on the wire 
surface will drift with velocity v towards the tip, and will be absorbed by it in an average 
time t ~ L/4v, where L is the wire length. Due to that, the wire radius will decrease 
roughly by 5R = — Xjl-nRLp. The free energy will first of all decrease by —l^jRp since 
an atom has left the wire for the (bulk-like) tip. However the wire radius decrease implies 
an additional free energy cost, all other atoms being left behind in a now thinner wire, 
amounting to — (2 / -f/R 2 p)(irR 2 Lp — 1)5R. This cancels one half of the previous free energy 
decrease (neglecting terms), so that the total free energy change for an atom leaving the 
wire and going to the tips is, to lowest order in 

AE = --y/Rp (2) 

Now the average velocity of surface atom, drifting from the wire middle towards the closest 
tip, approximately a distance L/4 away, will be v = p s (AE/ Ax), or v = 4rffi 3 / RLp. Assum- 
ing all surface particles within an atomic diameter 2r , whose number is N s = 47iLpr (R—r ), 
to drift in this manner, the radius will decrease with a thinning rate: 

AR = 6RN. = w a rj(l-r Q /R) 

At L/4v 3RL 2 1 ' 

where p = 3/4nrQ has been used. The predicted wire thinning rate calculated from this 
formula for wires of fixed length and variable radius, with a surface free energy ~ 
1.3 J /m 2 and diffusion constant D s ranging from 10~ 5 (typical of liquid Au) to 10~ 7 cm 2 /s, 
is shown in Fig. 1. As can be seen a typical value for the thinning rate can be in the order 
of 1 A/s for a nanometer thin wire, which seems of the correct magnitude to explain the 



experimentally observed thinning at room temperature ||20|| . The thinner and shorter a wire 
is, the faster it will further thin down, until it gets monatomic. This behavior should be 
amenable to future experimental verification. 

Once the thinning has reduced it to monatomic, the wire can only stretch under the tip 
suction, until it will eventually break. It is therefore mandatory to estimate the stability of 
a wire under stretch. 
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II. THE INFINITE MONATOMIC GOLD WIRE 



The electronic structure of an infinite, regular gold monatomic wire has been calculated 
as a function of the interatomic spacing. The calculation was performed within density 
functional theory in the LDA, and also in GGA, using the nonlocal ultrasoft pseudopotentials 
described in reference P5| . A repeated cell geometry was used, with a large inter- wire 
distance of 10.6 A, and a plane wave expansion with a cutoff energy of 25 Ry for the 
wavefunctions and 200 Ry for the charge density. The calculated total energy curves versus 
distance are shown in Fig. ^. They have minima of about 2.2 eV/atom (LDA), or 1.7 
eV/atom (GGA) |2|] at an equilibrium distance of about 2.50 A (LDA) or 2.55 A (GGA). 
For reference, with the same pseudopotentials and cut-offs, the calculated cohesive energy 
for bulk Au is 4.4 eV (LDA), 3.3 eV (GGA), against an experimental value of 3.9 eV. 
For comparison we also present in Fig. |2] the wire cohesive energy calculated using an 
alternative local orbital based LDA and GGA scheme, the SIESTA program [25,26], as well 
as those given by the empirical glue potential for gold P7|[ . The wire cohesive energy per 
atom is 2.8 eV (LDA), 2.1 eV (GGA), |2§ and only 0.4 eV (glue). The corresponding bulk 
cohesive energies provided by the local orbital calculation is 5.4 eV (LDA) 3.9 eV (GGA), 
3.9 eV (glue). The four ab-initio calculations are in substantial agreement (we have not 
attempted a further improvement by extending the basis) while the glue model is seen to 
severely underbind, and should be discarded for such an ultra-thin wire. |29[| Stretching the 
wire interatomic distance from its equilibrium value to, say, 3.5 A would cost, according 
to the local orbital GGA calculation, a stretching energy of about 4.0 ±0.5 eV, nearly 
exactly compensating the energy gain corresponing to one gold atom sucked by one of the 
tips, namely the bulk cohesive energy, which is 3.9 eV. This would seem a very reasonable 
balance, supporting the concept that tip suction is indeed responsible for the observed wire 
stretching. 

However, this reasoning is not further supported by the electronic structure calculations 
where we find that even a relatively modest stretching endangers the wire stability. Using 
for simplicity a four atom cell with periodic boundary conditions, we tested for longitudinal 
stability by calculating total energy changes when one interatomic spacing is extended from 
a to a', while forcing the three remaining interatomic distances to shrink so as to keep a 
total lenght of 4a (see Fig. ^). The energy increase obtained by increasing a' relative to 
a represents the energy barrier which prevents the wire from longitudinal breaking. To 
our surprise, we found (Fig. 3) that the barrier is positive at the equilibrium spacing of 
2.55 A and just above, but disappears very quickly when the wire is stretched, making it 



unstable below 3 A. [|30| This result is puzzling, in view of the exceptional stability observed 
experimentally. However, there is still a distinct possibility that the real 4-atom wire hanging 
between well-tapered tips could in many ways behave differently from the tip-less wire just 
studied. We thus decided to investigate further a more realistic case. 



III. SHORT MONATOMIC WIRE BETWEEN (110) TIPS: AB-INITIO 

SIMULATION 

For the first-principles MD relaxation, we generated a starting configuration in the fol- 
lowing manner. We started with an FCC bulk-ordered cubic-shaped piece of 240 atoms. 
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Keeping both extreme layers rigid, we separated them at constant velocity along the (110) 
direction and simulated the evolution of the neck using glue model MD (so long as atomic 
coordination not too poor, the glue model provides quite reasonable interatomic forces for 
Au.) The stretching velocity was ~ 0.14 mA/fs, and the temperature was kept very low. A 
neck formed and broke. After neck breaking we stopped the simulation and inserted three 
additional atoms in a row, forming a wire bridging the opened gap. Finally, we relaxed the 
new configuration, and the interatomic spacing of the relaxed wire reached about 2.9 A. We 
further pulled the tips, stretching the monatomic wire up to the very moment before it would 
break and stopped there, finally relaxing and checking for stability with glue forces. The 33 
central atoms of this last configuration, shown in Fig. f| form our initial configuration, with 
two tapered tips and a 4-atom wire with an interatomic spacing of about 3.5 A . 

Starting with this configuration, we calculated the electronic structure and the forces on 
all 33 atoms using the local basis LDA j3l| . With these forces we conducted a short ab initio 
MD relaxation of 40 steps of 5 fs each plus 60 steps of 7 fs each. The resulting ab-initio 
evolution after 100, 200, 340, 480 and 620 fs is also shown in Fig. |j. We see that interatomic 
distances in the wire change rapidly. However, while most of them shorten, especially close 
to the tip, one distance is growing with time, and the wire is clearly breaking. It appears 
therefore that within this approximation, the short wire hanging from the tips is unstable, 
very much as the infinite wire was. 



IV. DISCUSSION 

We believe we have understood the thinning and stretching of wires between tips, and 
provided a thinning rate dependence upon the wire length and radius that can be checked 
experimentally. The exceptional stability of the stretched monatomic wire remains instead 
mysterious since state-of-the-art calculations predict that it should break by far. While 
independent measurements would clearly be desirable in order to confirm the large spacings 



observed in reference [2D], we can presently only speculate on the possible reasons why our 
theory might be in error on this problem. 

A first possibility is that electron correlations, weak in bulk Au, could become strong 
in severely stretched wires, owing to the exceedingly poor coordination, and to the large 
interparticle distance. In particular, typical one dimensional instabilities, involving anti- 
ferromagnetism or dimerization could play a role. Of them, dimerization is allowed in our 
simulation, but magnetism is not. While this is clearly an avenue to be explored, it seems 
at this stage unlikely that it could alter the energetics very drastically. All the same, if the 
stretched wire turned magnetic (and presumably also insulating), the corresponding small 
energy decrease might be sufficient to generate a small change of curvature and therefore a 
small barrier against breaking. 

A second factor which we neglected throughout is long-range correlations, and partic- 
ularly van der Waals forces. Van der Waals attraction is known to be strong especially 



between blunt tips [J32|| . This attraction could be enhanced by the presence of a conduct- 



ing wire joining the two tips, thus providing the missing barrier. This effect is absent in 
LDA/GGA which does not include long-range dispersion forces. 

An alternative possibility is that the TEM electron beam could be influencing, by heating 
or charging, the wire, and causing its anomalous behavior. However both tips are electrically 



5 



grounded, and it is further estimated |33| that only about one electron per second collides 
elastically with the wire, the inelastic cross section being even smaller. 

Clearly the problem remains open, and more experimental and theoretical work will be 
needed to solve it. 

We are grateful to Prof. K. Takayanagi and his group for enlightening discussions. 
We acknowledge financial support from INFM/G (Project " Nanowires" ) , INFM/F (PRA 
LOTUS), by MURST through COFIN97, and Spain's DGES, project PB95-0202. Work by 
J.A.T. was carried out under TMR grant ERBFMBICT972563. 
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FIGURES 



FIG. 1. Thinning rate of tip suspended 20 A long gold wires, from Eq. 2, calculated for 
different surface diffusion constants. 



FIG. 2. Cohesive energy per atom for an infinite monatomic Au wire. 



FIG. 3. Energy variation upon attempted wire breaking. The wire is not stable when the 
interatomic spacing is stretched above 2.8 A. 



FIG. 4. Snapshots of the ab-initio MD relaxation for increasing time. The spontaneous 
breaking of the wire is evident. 
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